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Abstract. 

We study luminosities of millisecond pulsars in globular clusters by fitting the observed lumi- 
nosity distribution with single and double power laws. We use simulations to model the observed 
distribution as the brighter part of some parent distribution for Terzan 5 and try to find a model 
which simultaneously agrees with the observed diffuse radio flux, total predicted number of pul- 
sars and observed luminosity distribution. We find that wide ranges of parameters for log-normal 
and power-law distributions give such good models. No clear difference between the luminosity 
distributions of millisecond pulsars in globular clusters and normal disk pulsars was seen. 
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INTRODUCTION 

Millisecond pulsars (MSPs) in globular clusters (GCs) can be used as a tool to under- 
stand the properties of GCs as well as the recycling process in the dense stellar environ- 
ments inside GCs [1]. As luminosity is a fundamental property of pulsars, one necessary 
step to achieve this goal is to understand the luminosity distribution of MSPs in GCs. A 
full dynamical approach, where one models evolution of pulsars and observational lim- 
its following appropriate choice of birth distributions of pulsar parameters, can not be 
adopted for GC pulsars where it is difficult to model the effects of stellar encounters and 
the cluster potential. The simplest way is to use a snapshot approach where one mod- 
els pulsar luminosities as observed. We first adopt this method and then, using Monte 
Carlo simulations, try to find a good model which not only fits the observed luminosity 
distribution, but also agrees with the total radio and y-ray fluxes for Terzan 5. 

ANALYSIS AND RESULTS 



We consider GC pulsars with known flux values and spin periods smaller than 100 ms 
and exclude the GCs having less than four such pulsars. Thus our dataset contains 79 
pulsars in 9 GCs. We use the flux values measured at 1400 MHz (5i4oo) if available, 
otherwise we estimate the value of 5i4oo setting the value of the spectral index (a) to 
be —1.9, which is the mean value obtained from the MSPs in GCs having flux values 
reported at multiple frequencies. This value agrees with the earlier estimate for 19 
MSPs including two GC MSPs [2]. "Pseudo-luminosities" (L1400) of the pulsars have 
been calculated using the relation L1400 = Sx^qqcI^ where d is the distance of the host 
GCs from the sun. We confirm the earlier conclusion [3] that the choice of a in a 



realistic range does not significantly affect the complementary cumulative distribution 
(CCD) of pulsar luminosities. With these data, we first revisit the study of Hui et al. 
(2010; hereafter HCTIO) [4] who modeled CCDs of luminosities of GC pulsars as 
A^(> ^i40o) = A^o^uoo- We find A^o = 59± 1, ^ = -0.80±0.03 when we fit all 79 pulsars 
and A'^o = 74^4, q = —1.06 ±0.06 when we fit 48 pulsars with L1400 > 1.5 mJy kpc^; 
whereas HCTIO found A^o = 68 ±2, ^ = -0.58 ±0.03 when they fitted all 78 pulsars and 
A^o = 91 ±6, = -0.83 ±0.05 when they fitted 58 pulsars with L1400 > 1-5 mJy kpc^. 
The fit improves if we keep only pulsars with L1400 > 1.5 mJy kpc^. Moreover, we find 
that the fit significantly improves if we adopt two power laws. Fig. 1 shows the single 
and double power law fit of GC pulsars with L1400 > 1.5 mJy kpc^. We also fit a single 
power law for pulsars with L1400 > 0.5 mJy kpc^ in individual GCs separately and find 
that our fitting parameters are different from those obtained by HCTIO shown in Table 
1. This is due to two reasons. Firstly, our dataset is somewhat different. HCTIO did 
not exclude pulsars with spin period greater than 100 ms, they kept the pulsars in NGC 
6441 which we exclude, we have 7 pulsars in M15 but HCTIO did not consider this GC 
at all. Secondly, we use recent distance estimates of GCs which are usually different 
from those used by HCTIO, the largest difference is for Terzan 5. 
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FIGURE 1. Single (left) and double power-law (right) fit of GC pulsars with L1400 > 1-5 mJykpc^. 
Double power-law fit parameters as Nqj ~ 70^g, qi = —0.97 ±0.13 for L1400 < 4.0 mJy kpc^ and A^o./i = 
134;^?, qi, = -1 .40 ± 0.21 for L1400 > 4.0 mJy kpc^. 

Similarly, one can use any other distribution function to fit the observed CCD, but 
instead of doing so, we invoke a better method which simultaneously agrees with the 
observed diffuse radio flux {So,tot), total predicted number of pulsars (Ntot) and the 
observed luminosity distribution. The luminosities of Ntot pulsars are simulated from 
a chosen distribution and the distribution of luminosities of A^^ pulsars having L > Lm,o 
is compared with the observed distribution using KS tests where the KS probability i^^. 
should be high for a good model. Lm,o is the observed minimum luminosity. But there 
can be models which give large/small values of Ns in comparison to the observed number 
of pulsars {No), but the CCD having the same shape as the observed CCD. This can not 
be a good model although will be large. To overcome this problem, we check whether 



TABLE 1. Power Law fit parameters for different GCs 
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the "goodness factor" X = is ~ 1. Here we define 

^ iSs,tot — So,tot)^ . (Ns—No)'^ U C V^' C /1\ 

Q = — 7^ ^ + - — , wheTeSs,tot^2^Si (1) 

{^o,tot) i=0 

for simulated fluxes 5;. Presently we apply this method only in the case of Terzan 5 for 
which Ntot has been predicted to be ~ 1 80 from the total y-ray flux [5] and So,tot was 
estimated to be 5.22 mJy [6]. We use log-normal, power-law and exponential distribution 
functions over a wide range of parameters. We plot P^s and X for all three distributions 
in Fig. 2. We plot the parameters of log-normal and power-law distributions along x and 
y axes and the values of Pks and X as color codes. For the single-parameter exponential 
distribution, we show the value of the parameter (mean~^) along the x axis and the 
value of and X along the y axis. For log-normal and power-law distributions, we 
find wide ranges of parameters where both P^s and X have high values, so it is difficult 
to pinpoint good luminosity distribution functions. But for the exponential distribution 
we do not see a parameter space where both and X have high values and conclude 
that the exponential distribution can not describe pulsar luminosities adequately. In 
summary, we find that both log-normal and power-law distributions can reproduce the 
observed luminosities of MSPs in Terzan 5 but it is difficult to find exact parameters for 
a good model. Moreover, the model invoked [7] for normal disk pulsars (FK06 model, 
log-normal distribution with mean -1.1 and standard deviation 0.9) falls in the good 
region (Fig. 2), suggesting that GC MSPs might have the same luminosity distribution 
as normal disk pulsars. It will be interesting to perform similar studies for other GCs 
which may constrain the distribution functions leading to the answer of the question 
whether GC pulsars form a different population from disk pulsars. 
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FIGURE 2. Plot of ffa and X for log-normal, power-law and exponential distributions. For log-normal 
and power-law distributions, we plot the parameters along the x and y axes and the value of Pj^s or X 
as color codes. The parameters for the log-normal distribution are mean (jj.) and standard deviation (cj), 
and the parameters for the power-law distribution are minimum luminosity (L„,) and index (j3). For the 
exponential distribution, we show the value of the parameter (A) along the x axis and the value of or 
X along the y axis. 
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